Many-spin effects and tunneling splittings in Mni 2 magnetic molecules 



O 

o 

(N 
>v 



(N 



H. A. De Raedt and A. H. Hams 1 and V. V. Dobrovitski, M. Al-Saqer, M. I. Katsnelson, B. N. Harmon 2 

1 Institute for Theoretical Physics and Materials Science Centre, 
University of Groningen, Nijenborgh 4, NL-9747 AG Groningen, The Netherlands 
2 Ames Laboratory, Iowa State University, Ames, I A 50011, USA 

We calculate the tunneling splittings in a Mni2 magnetic molecule taking into account its internal 
many-spin structure. We discuss the precision and reliability of these calculations and show that 
restricting the basis (limiting the number of excitations taken into account) may lead to significant 
error (orders of magnitude) in the resulting tunneling splittings for the lowest energy levels, so that 
an intuitive picture of different decoupled energy scales does not hold in this case. Possible routes 
for further development of the many-spin model of Mni2 are discussed. 

PACS numbers: PACS: 75.40.Mg, 75.50.Xx, 75.10.Dg, 75.45.+j 
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I. INTRODUCTION 

Molecular mag netsH have proven to be very suitable 
systems for the study of mesoscopic tunneling effects 
in magnetic materials. A number of impressive ex- 
perimental results have been ghtained recently, such 
as thermally-assisted tunnelingatj, ground state - to - 
ground state tumaelingaia and topological phase effects 
in spin tunnelingp. Among others, the molecular mag- 
net Mn 12 Oi2(CH 3 COO) 16 (H20)4 (below referred to as 
Mni2) has received special attention. The effect of res- 
onant magnetization tunneling has n been first observed 
and studied in detailed experimentaj'u on Mni2 , and, at 
present, a substantial amount of reliable experimental 
data has been collected. Quantitative analysis of these 
experiments is a challenging theoretical problem involv- 
ing fundamental issues about tunneling phenomena in 
mesoscopic magnetic systems. The basic prerequisite for 
solving this problem is our ability to evaluate accurately 
and reliably the energy splittings occuring as a Result 
of tunneling between two (quasi)degenerate levelal At 
present, carefully designed magnetic relaxation experi- 
ments at low and ultralow temperatures (tens or hun- 
dreds of milliKelvins) can detectnm the changes in relax- 
ation time caused bv-,the splittings of order 10 _2 -10 -4 
K, and even smallerQ, of order 10~ 6 -10~ 7 K. The re- 
laxation time data obtained in these experiments give 
information (although indirect) about the splitting val- 
ues, so that predictions of the theoretical models can be 
compared with experimental results. 

Conventionally, the molecular magnet Mni2 is consid- 
ered as a large single spin S = 10 with quasidegenerate 
levels S z = +M and S z = —M split because of tunneling. 
However, the single-spin Hamiltonian is a phenomeno- 
logical construct; in reality, this is a many-spin system, 
consisting of 12 manganese ions coupled by exchange in- 
teractions. Here, using Mn 12 as a well-studied example, 
we address the problem of reliable many-spin calculation 
of the tunneling splittings in molecular magnets. Such 
a calculation is a very complicated task. For example, 
the Hilbert space of the spin Hamiltonian describing a 
molecule of Mni2 consists of 10 s levels, while the small- 



est tunneling splittings in Mni2 are of order of 10~ 10 
Kelvin (as measured in RefJa for m — ±10). The brute- 
force direct calculation of tiny tunneling splittings in this 
system, even for several low-lying states, is beyond the 
capabilities of modern computers. The general strategy 
to solve this problem is to truncate the full Hilbert space 
thus reducing consideration to a much smaller number of 
relevant energy levels. This idea, implemented in a rather 
sophisticated way, forms a basis of several approaches for 
the evaluation of tunneling phenomena, such as quantum 
Monte-Carlo methodsu, stochastic diagonalizationli3, and 
instanton calculationsO. 



To our knowledge, all calculations of the tunneling 
splittings in molecular magnets starting from realistic 
Hamiltonians have employed truncation of the Hilbert 
space in a much more straightforward, and much less 
justified manner. High-energy basis states, assumed to 
be irrelevant, are being explicitely excluded from consid- 
eration, and only the low-energy part of the spectrum 
is being taken into account! 1 ! In the present paper, we 
calculate tunneling splittings using the many-spin model 
of Mni2, examining the accuracy and reliability of this 
straightforward scheme. We demonstrate that, because 
of strong Dzyaloshinsky-Morya interactions present in 
Mnl2, the splitting values obtained in this way are unre- 
liable. We also consider the sensitivity of the calculated 
splitting values to variation in the Hamiltonian parame- 
ters, and determine the accuracy needed for reliable split- 
tings calculation. 



The paper is organized as follows. In Section II, we 
discuss the 8-spin model of Mni2 and the methods used 
to calculate tunneling splittings based on this model. We 
also consider the stability of the results with respect to 
possible limitations of the model Hamiltonians. In Sec- 
tion HI we consider the reasons for the failure of the 



energy-based truncation scheme in the splittings calcula- 
tions. Our conclusions can be found in the Summary. 



2 



II. 8-SPIN MODEL OF Mn 12 AND 
CALCULATIONS OF THE TUNNELING 
SPLITTINGS 

The cluster M1112 consists of eight Mn 3+ ions having 
spin 2 and four Mn 4+ ions having spin 3/2, coupled by 
exchange interactions. The total number of spin states in 
Mni2 is 10 8 , and a corresponding Hamiltonian matrix is 
rather large to be treated by modern computers. To over- 
come this difficulty, we can employ the natural hierarchy 
of interactions present in Mni2. The antifcrromagnetic 
exchange interactions J\ ~ 220 K between Mn 3+ axid 
Mn 4+ ions are significantly stronger than all the othersEj, 
so corresponding pairs of Mn 3+ and Mn 4+ ions can be 
considered as stiff dimers with the total spin s = 1/2, 
thus giving rise to the 8-spin model of Mni2 . The range 
of validity of the 8-spin model, and the corresponding 8- 
spin Hamiltonian of Mn 12 have been considered in Ref. 
[Ti] . After examination of different possible interactions, 
the following Hamiltonian has been proposed: 

i (k,l) i=l 

(i,j) 

Here, Si and are the spin operators for the large spins 
S = 2 and small dimer spins s = 1/2, correspond- 
ingly (the subscript i indexes the spins). The first two 
terms describe isotropic Heisenberg exchange between 
the spins. The third term describes the single-ion easy- 
axis anisotropy of large spins. The fourth term represents 
the antisymmetric Dzyaloshinsky-Morya (DM) interac- 
tions in Mni2, where D ZJ is the Dzyaloshinsky-Morya 
vector describing the DM-interaction between i-th small 
spin and j-th large spin. Existence_pf DM-interactions 
in Mni2 has been suggested in Ref.E3, and their magni- 
tude has been estimated in Ref.EfJ based on the neutron 
scattering datalia. The molecules of Mni2 possess a four- 
fold rotational-reflection axis (symmetry S4) imposing re- 
strictions on the DM- vectors D 4J , so that Dzyaloshinsky- 
Morya interactions can be described by only three param- 
eters D x = Dl' s , D y = #i<*Land D z = D 1 ^. 

It has been demonstrated^ that the above model sat- 
isfactorily describes a rather wide range of experimen- 
tal data, such as the splitting of the neutron scattering 
peaks, results of EPR measurements and the temperature 
dependence of magnetic susceptibility. Here,for calcula- 
tions we use the parameter set A from Ref. [lif 



set A: 

J = 0, J' = 105 K, K z = 5.69 K 
D x = 25 K, D y = 0, D z = -1.2 K. 



(2) 



molecule). However, this set of parameters should not 
be considered as being accurately determined, since the 
amount of the experimental information available is not 
yet sufficient to achieve particularly reliable parameters. 
In the Hamiltonian (Q), only the fourth term, repre- 
senting the Dzyaloshinsky-Morya (DM) interactions, can 
lead to tunnelingli3: the first three terms conserve the z- 
projection of the total spin S z and can not induce tunnel- 
ing between levels with different S z , while the DM-term 
mixes levels with different S z . In what follows, we will 
label the energy levels by the value of S z . Although it is 
not an exact quantum number, we can formally consider 
the DM-interaction as a perturbation, and use perturba- 
tion theory terminology. 

The following values of the tunneling splittings corre- 
sponding to the parameter set (||) have been obtained by 
the diagonalization of the full Hamiltonian matrix (of the 
size 10 4 x 10 4 ) using quadruple precision arithmetics: 



AS(±10) 
AE{±6) 
AE(±2) 



1.18- 10" 15 K, 
3.87- 10~ 8 K, 
4.17- 10" 2 K. 



AE{±8) 
AE(±A) = 



= 1.06- 1CT U | 
2.08 • 10~ 6 K, 



which also gives a good description of the response of 
Mni2 molecules to a transverse magnetic field (exter- 
nal field applied perpendicular to the easy axis of the 



The splittings for odd values of S z are not shown: they 
constantly remain at the level of the numerical precision 
of the calculations (of order of 10~ 19 K)E3. In Mni2, 
these splittings should be zero since the fourfold symme- 
try of the molecule imposes certain restrictions on the 
symmetry of the spin Hamiltonian and makes some ma- 
trix elements vanish. In the single-spin model of Mni2 
this property of the spin Hamiltonian is introduced ex- 
plicitly, by retaining only those operators which possess 
the required fourfold symmetry. In the many-spin simu- 
lations, we obtain the same result independently. 

The first question to pose concerns the accuracy of 
the level splitting evaluation. Parameters of the Hamil- 
tonian are determined with some finite precision, and a 
small error (say, of the order of several Kelvin) affects 
the level energy by an amount of order of Kelvin, which 
is much larger than the very small value of tunneling 
splitting (of order of 10 -12 K). Does it deprive the calcu- 
lational results of all meaning? To answer this question, 
we note that the levels \S Z = +AI) and \S Z = —M) are 
degenerate due to exact symmetry properties of the spin 
Hamiltonian, and, in the absence of the DM-term, would 
be degenerate at any value of parameters. Therefore, 
the tunneling splittings AE + m.-m are governed only by 
the strength of the interaction which breaks the symme- 
try, i.e. the DM-interaction. If the parameters of the 
Hamiltonian are determined with reasonably small rela- 
tive error, and if the numerical calculation is done with 
sufficient precision, then the relative error of the level 
splittings will also be small. This conclusion is supported 
by our calculations: a 10% variation in the Hamiltonian 
parameters leads to the variation in the splitting values at 
most by a factor of ten, so that accurate determination of 
the Hamiltonian parameters is necessary for reliable cal- 
culation of the tunneling splittings. If only a logarithmic 
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accuracy in the splitting values is needed, then the 10% 
uncertainty in the Hamiltonian parameters is sufficient. 

However, there is another, much more important 
source of possible error. The description of the Mni2 
molecule by the 8-spin model requires a full, high- 
precision diagonalization of the Hamiltonian matrix with 
dimensions 10 4 x 10 4 . Solving this problem is rather time- 
consuming. Matrices of that size can be processed very 
effectively using Lanczos-type methods, but the applica- 
tion of these methods to the tunneling splitting calcu- 
lations constitutes quite a difficult problem by itself. A 
very large number of iterations is needed to achieve the 
necessary precision and in addition the precision is hard 
to control when the level separation is very small, so that 
special techniques are necessary. 

Therefore it is natural first to explore another ap- 
proach, namely, to omit high-energy basis states, retain- 
ing only the low-lying part of the spectrum where basis 
levels have energies less than some threshold value E cut . 
This approach has been adopted extensively and in fact, 
we are not aware of any calculations of tunnel-splitting 
of magnetic molecules done in a different way: calcula- 
tions based on both the single-spin and the many-spin 
modety have employed this method. In this paper, we 
assess the validity of this energy-based truncation ap- 
proach by considering the dependence of the tunneling 
splittings AE + m,-m f° r different pairs of degenerate lev- 
els \S Z = +M) and \S Z = —M) on the number of low- 
est levels Mow actually used in calculations (or, in other 
words, their dependence on the energy threshold E cnt ). 

A brief description of the basis states is in order. We 
first consider the first two exchange terms in the Hamil- 
tonian of Eq. (Q) and diagonalize within the manifold of 
all the 8-spin configurations yielding states with S z — 0; 
there are 1286 energy eigenvalues corresponding to eigen- 
vectors with S ranging from to 10. The distribution of 
states is: (10,1), (9,7), (8,24), (7,56), (6,104), (5,164), 
(4,220), (3,248), (2,232), (1,168), (0,62), where the first 
number in parenthesis is the value of S and the second 
is the number of levels with this value of S. With the 
25 + 1 degeneracies included, there are exactly 10000 
states. These are the basis states which are then used 
to diagonalize the full Hamiltonian, including anisotropy 
and DM terms. 

The initial increase in the number of basis states con- 
sidered, 7Vi ow , leads to an overall increase in AE + m,-m 
accompanied by oscillations (see Fig. Q). After N\ ov/ 
achieves the value of about 700, the oscillations have 
become small and AE + m,-m versus N\ ow exhibits a 
plateau. This saturation lead in Ref. [ll] to the conclu- 
sion that the resulting values give the actual splittings 
with sufficient accuracy. But this conclusion is wrong. A 
further increase of the number of levels leads to a resur- 
rection of the oscillations at N\ ow ~ 1200, with a quite 
pronounced jump in AE + m.-ai for Ai ow ~ 1700. For a 
larger number of levels, the situation repeats itself: the 
values of the splittings reach another plateau, then oscil- 
lations appear again with a subsequent jump, etc. We 



have traced this behavior up to N\ ow ~ 3000, which is 
already 1/3 of the total number of levels. The observed 
behavior of AE + m,-m is, in our opinion, a very clear 
signal that energy-based truncation of the Hilbert space 
is not a good strategy for the computation of tunneling 
splittings: it gives unreliable results. 

The rather sharp jumps in the tunneling splittings as 
discussed above and illustrated in Fig. [l] are associated 
with the inclusion of basis states with large S values. Be- 
cause of the selection rule for the DM term (S — > S ± 1), 
the >S = 10 ground state only couples with S = 9 states. 
States with smaller S values affect the splittings more 
indirectly by coupling with other states which eventually 
couple to the ground state. While the states with large 
S cause jumps in the splitting values, there are few of 
them, and the smaller coupling of smaller S states still 
is significant because of the cumulative effect of so many 
states (see the distribution given above). Therefore, the 
evaluation of tunneling splittings for a general system 
possessing strong DM interactions requires consideration 
of sufficiently large portion of Hilbert space. 

It is noteworthy that the same truncation method 
works rather well for calculations of the energies of well- 
separated levels. To compare the model against most 
of the experiments, it suffices to know the positions of 
the levels with much less precision, usually an error less 
than 0.1 K is already adequate. This level of precision 
can be obtained by taking into account A^ ow ~ 1000 lev- 
els (i.e., 1/10 of the total Hilbert space). Even using 
N\ ow ~ 500, the error in the level position is less than 1 
K even for the states of energy about 60 K. Therefore, 
the matrix-truncation approach is adequate for fitting the 
model parameters to experimental data. But the calcu- 
lations of the tunneling splittings should be done using 
the full Hamiltonian matrix. 



III. DISCUSSION 

We have shown that truncating the Hilbert space leads 
to large errors in the calculated values of tunneling split- 
tings. But actually, any sensible Hamiltonian is in- 
evitably obtained due to some truncation of the Hilbert 
space. For example, the Hamiltonian ([!]) epi be con- 
sidered as a result of the two-step proceduretS (i) pro- 
jection of the real many-electron Hamiltonian onto the 
subspace of suitably chosen single-electron orbital states, 
yielding a general spin Hamiltonian of the molecule; and 
(ii) projection of the resulting spin Hamiltonian onto the 
subspace of the 8-spin model. This procedure is usu- 
ally justified (at least, at the heuristic level) by invoking 
some kind of perturbation or WKB-theory arguments, 
and corresponds to an intuitive idea of different, practi- 
cally independent energy scales. 

However, in the case of the tunneling splittings, we 
see that very different energy scales significantly affect 
each other. Why do the same arguments not work if we 
truncate the 8-spin Hamiltonian? In our opinion, this 
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takes place because the conditions of the applicability of 
WKB-reasoning (or similar arguments based on pertur- 
bation theory) are not satisfied. The spin of the system 
S = 10 is too small, so that the instanton actional on 
the trajectories corresponding to the 8-spin model is not 
large enough. Indeed, for systems with well-separated 
levels, the quasiclassical approximation usually already 
works reasonably for a-total spin S ~ 2-3. However, as 
has been demonstrated!!^, to apply the same type of ar- 
guments to the splitting calculations, the (normalized) 
instanton action Sj should exceed the value of 12. For 
the model employed in Ref. |l^, this corresponds to the 
system with a total spin (more exactly, with the total 
antiferromagnetic vector) of order of several thousand. 
Thus, the tunneling splittings, in general, appear to be 
much more sensitive to the method of calculation than 
the level energies themselves, and conditions for applica- 
bility of the conventional WKB-reasoning are consider- 
ably more stringent (though for Mni2 they can of course 
be different from the condition Si > 12). Qualitatively 
this agrees with our observations (see Section II) . Even a 
rather severe truncation of the Hilbert space has a minor 
effect on the level energies, while correct values of the 
tunneling splittings require a diagonalization of the full 
Hamiltonian. 

Briefly, these arguments can be expressed in a rather 
obvious form: the 8-spin model is not "macroscopic 
enough" to justify the truncation of the Hilbert space 
by some WKB or similar perturbation approach. In this 
case the intuitive picture of different independent energy 
scales is misleading. 

This conclusion raises important questions, namely, is 
the 8-spin model, being the result of the truncation of, 
e.g., 12-spin Hamiltonian, sufficient to predict reliably 
the tunneling splittings (or, in other words, is the 12-spin 
model "macroscopic enough" to be truncated)? What is 
the minimal model allowing the splittings to be calcu- 
lated correctly? We believe that these are key questions, 



not only for Mni2 but for the whole class of magnetic 
molecules. For this purpose, ab-initio calculations of the 
exchange and anisotropic intramolecular interactions in 
Mni2 could be very useful. Also, reliable experimental 
data for the tunneling splittings would obviously be of 
great value for further development. 



IV. SUMMARY 

We have calculated the tunneling splittings ha— Mni2 
on the basis of the 8-spin model proposed earlieio. We 
have shown that rather accurate knowledge of the Hamil- 
tonian parameters is needed for the accurate splitting 
calculations; although, for logarithmic accuracy, 10% er- 
ror in the parameters can be tolerated. Furthermore, 
we have demonstrated that a reliable calculation of the 
tunneling splittings for a system with strong DM inter- 
actions requires the use of the full Hamiltonian matrix. 
We have explicitely shown that an energy-based Hilbert 
space truncation scheme can be successfully used for the 
determination of the level energies, but leads to erroneous 
results when applied to the splitting calculations. 
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